Beyond power laws: Universality in the average avalanche shape 
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We report the measurement of multivariable scaling functions for the temporal 
average shape of Barkhausen noise avalanches, and show that they are consistent 
with the predictions of simple mean-field theories. We bypass the confounding 
factors of time-retarded interactions (eddy currents) by measuring thin permal- 
loy films, and bypass thresholding eff"ects and amplifier distortions by applying 
Wiener deconvolution. We find experimental shapes that are approximately 
symmetric, and track the evolution of the scaling function. We solve a mean- 
field theory for the magnetization dynamics and calculate the form of the scaling 
function in the presence of a demagnetizing field and a finite field ramp-rate, 
yielding quantitative agreement with the experiment. 

The study of critical phenomena and universal power laws has been one of the central ad- 
vances in statistical mechanics of the second half of the last century, in explaining traditional 
thermodynamic critical points avalanche behavior near depinning transitions , and 
a wide variety of other phenomena Q]. Scaling, universality, and the renormalization group 
claim to predict all behavior at long length and time scales asymptotically close to critical 
points. In most cases, the comparison between theory and experiments has been limited 
to the evaluation of the critical exponents of the power law distributions predicted at crit- 
icality. An excellent playground of scaling phenomena is provided by systems exhibiting 
crackling noise, such as the Barkhausen effect in ferromagnetic materials |5|. Here we focus 
on the average functional form of the noise emitted by avalanches — the temporal average 
avalanche shape {4]. 

This avalanche shape has been measured for earthquakes Gj and for dislocation avalanches 
in plastically deformed metals jl, 8|, but the primary experimental and theoretical focus has 
always been Barkhausen avalanches in magnetic systems j^,y,[9-ll|. Theory and experiment 
agreed well for avalanche sizes and durations, but the strikingly asymmetric shapes found 
experimentally in ribbons disagreed sharply with the theoretical predictions, for which 



the asymmetry in the scaling shapes under time reversal was at most very small |6[. 
(We note that the relevant models are not microscopically time-reversal invariant; temporal 
symmetry is thus emergent). Doubts about universality [4] were resolved when eddy currents 
were shown to be responsible for the asymmetry, at least on short time scales 12j, but the 
exact form of the asymptotic universal scaling function of the Barkhausen avalanche shape 
still remained elusive. 
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In this paper, we report an experimental study of Barkhausen noise in permalloy thin 
films, where a careful study of the average avalanche shapes leads to symmetric shapes, 
undistorted by eddy currents (which are suppressed by the sample geometry). We provide 
a quantitative explanation of the experimental results by solving exactly the mean-field 
theories for two general models of magnetic reversal: a domain wall dynamics model |[13i] and 
the random-field Ising model [l4|. The two mean-field theories are shown to be equivalent, 
and allow us to compute the average avalanche temporal shape as a function of typical 
experimental control parameters such as the field rate and the demagnetizing factor. 

The relevant statistical information encoded in the Barkhausen noise could in principle be 
derived from the joint two-point time- velocity distribution Gc^kiV, t; V, t + A), yielding the 
conditional probability that the noise at time A is equal to V if it was equal to V at time 
t. Here c is the external field rate and k the demagnetizing factor. The standard avalanche 
distributions can be derived from this parent distribution, e.g. the duration distribution is 
given by 

/•oo 

P(T) = - / dVd^G,^k{0^,t;V,t + A)\^=T. (1) 
Jo 

The renormalization group makes use of an emergent scale invariance for Barkhausen noise. 
Here, the two-point time-velocity distribution, when coarse-grained in time by a factor b 
and rescaled downward in velocity by a factor b^, rescales at long durations to itself: 

G,4V,t;V',t + A) =b-'^G,-,{V/b\t/b; V /b^,{t + A)/b), (2) 

where c = b^c is the rescaled field rate, k = h^k the rescaled demagnetizing factor and x, w 
are universal scaling exponents. Repeating n rescalings until A/6" = 1 leads to a universal 
scaling form 

a,fc(V^, t- V, t + A) = A-^^g{A-^V, A-W, A^'c/co, A^'fc/fco). (3) 

where ^ is a universal multivariable scaling function and cq and ko are the small scale values 
of the field rate and demagnetizing factor, respectively. 

Universal scale invariant forms can then be derived for all statistical quantities of interest 
including the temporal average shape. For avalanches of duration T, the average shape is 
defined as the average velocity for avalanches that begin and end at \^ = in a duration T. 
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It has an associated universal scaling form, 

(K,fc(A|T)) = J dV'V'G,,k{0-',0;V',A)G,.k{V',A;0,T)/ J dV'G,,k{0'' ,0;V' , A)G,,k{V' , A;0,T) 
= T^V{A/T,ik/ko)T^,c/coTy), (4) 

where V{X, K,C) is a universal scaling function, dependent on the rescaled time A = A/T, 
the rescaled demagnetizing factor K = {k/kQ)T^, and the driving field C = (c/co)T^. It is 
universal in the sense that in the scaling regime does not depend on microscopic features of 
the material, so long as the system is at a critical point. 

We record the Barkhausen noise by a standard inductive technique on a 1 /im thick 
permalloy film, with polycrystalline structure (see methods section for details on the mea- 
surement and on sample preparation). The Barkhausen noise is composed by a series of 
intermittent pulses, due to avalanches in the magnetization, combined with background in- 
strumental noise. In most crackling noise phenomena, avalanches are usually identified by 
setting a threshold above the background noise V^h- This method works well if the signal 



to noise ratio is high, but can induce spurious effects otherwise [15|. Our thin films have a 
correspondingly weak signal, making the Vth inappropriate. In addition, the measurement 
apparatus has a response function that distorts the original pulses. We instead extract the 
pulses using Wiener deconvolution 1^, which optimally filters the background noise and 
bypasses the use of thresholds (see Fig. ISTI) . 

The universality class of the Barkhausen noise in a sample is usually identified by mea- 



suring the voltage distribution, the power spectrum and the distributions of avalanche sizes 
S and durations T ^. The present experiments display characteristic features of the mean- 
field universality class: we observe a field-rate dependence of the voltage distribution and of 
the avalanche size and duration distributions that is in excellent agreement with mean-field 
predictions (see Fig. [S2|a), (b) and (d), noting in each case the three variable scaling im- 
plied by the field rate and demagnetizing factor). Furthermore, the power spectrum and the 
conditional average size of an avalanche of duration T are rate independent and consistent 
with mean-field predictions (see Fig. [S2](c)). Finally, we focus on the measurement of the 
average avalanche temporal shape, considering all the avalanches of a given duration T and 
averaging the voltage signal at each time step t. (In practice, we average the avalanches 
using duration bins centered at T and with geometrically increasing sizes). Fig. ISST a) shows 
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the resulting nearly symmetric temporal average avalanche shape, which starts out parabolic 
and then flattens as the duration of the avalanche increases (cf. Fig. ISST c)). 



It has been argued that dipolar magnetic flelds are sufliciently long-rangec 
fleld theory should be quantitatively applicable in three dimensional samples 



that mean- 



18|. Similar 



considerations apply to many of the systems exhibiting crackling noise, such as dislocation- 



mediated plasticity [19] and earthquakes 20|] where long-range interactions are provided by 
elastic strains. Our fllms are thinner than polycrystalline ribbons known to exhibit mean- 
fleld behavior j4|, but thicker than previously studied 2D fllms j^. The classic mean-fleld 



theory for domain wall depinning is the single-degree of freedom ABBM model [13|, which 
treats the domain wall as a rigid object at position X{t), advancing in a random pinning 
fleld statistically chosen as a random walk in X, 

= ct-kx + w{x), (5) 

with {W{X)) = and {W{X)W{X')) = \X - X'\ (Brownian noise) and where T is the 
damping coefficient. The ABBM model predicts that the avalanche size and duration dis- 
tributions decay as power laws with rate dependent exponents |lo|. The exact form of the 
scaling functions for these distributions, including the cutoff to the power law behavior, can 
be computed in the limit c = [l^. To compare with experiments for c > 0, we resort to 
numerical integration of Eq. 13 obtaining close agreement, as displayed in Fig. [S2l 

The average avalanche shape in the ABBM model has not been extensively explored 



numerically, but an approximate analytical calculation [10|, l22| gave a lobe of a sinusoid 
V(A) ~ sin(7rA). In that calculation, the approximation {V{X{t),T)) ~ V{{X{t)),T)) was 
made. (That is, the fluctuations in the growth of the avalanche size with time were neglected, 
leading to a slight distortion of the predicted temporal average avalanche shape.) We avoid 
this approximation by using a transformation due to Bertotti [23] to a time-dependent noise 
with a variance proportional to the avalanche velocity V 

^ = c-kV + Vwm, (6) 

which can be solved in the Stratonovich interpretation in the limit c = 0, k = 0. After 
defining a new variable Y = V^^'^, we have a stochastic equation that we may solve explicitly. 
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Utilizing the resulting probability for the first-return to the origin, we calculate the temporal 



shape 



24i\ (cf. Supplementary Information) : 



(K=o,.=o(t|T)) = (l - ^) , (7) 

where in the future we will use t instead of A to denote the time elapsed from the beginning 
of the avalanche. We have numerically verified that the average shape is indeed an inverted 
parabola (and not the lobe of a sine wave [lQ{]). 

Interestingly, an inverted parabola was reported in numerical studies of the (nu- 



cleated) random-field Ising model in the mean-field limit (the 'shell' model) |14|. Front 
depinning and nucleated transitions have different upper critical dimension and different 
short-range critical exponents, and it was conceivable that their mean-field theories had 
different average shapes, albeit sharing critical exponents. 

Upon closer examination, these two models are the same in the continuum limit. The 
shell model has a set of interacting spins Mi = ±1 with random fields hi distributed by a 
Gaussian, interacting with a strength J/N with all other spins; each spin flips when the net 
field it feels, hi + H(t) + J/N'^j Sj is positive. Here, H(t) = Hq + ct is the external field, 
increasing with rate c, and the spin feels the magnetization M = (l/N) 5*-,- both through 
the infinite-range coupling J/N and the demagnetizing factor k. A shell- model avalanche 
proceeds in parallel, with a shell of V„ unstable spins fiipping at the n— th time-step, then 
triggering a new set of spins Vn+i to fiip: 

K+i = P(2JK), (8) 

where P is the Poisson distribution for the set of spins Vn+i to be included in the range 
{/, / + 2JVn}. For large V, we may approximate the Poisson distribution as a Gaussian 
P{V) = V + VV^, leading to Vn+i — Vn — —kVn + V JVn^ (t) which is clearly a discretized 
version of Eq. [6l where = 1 — 2 J is the distance to the critical threshold. This equivalence, 
in retrospect, provides an alternative explanation for the origin of Brownian noise statistics 



for the ABBM domain wall potential 18 1. 



Finally, we consider the effects of the two main physical perturbations, the driving field 
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c and the demagnetizing factor k. The transformed Eq. |6]now takes the form 

Under rescahng \/V = Y b~^^'^Y, t — )■ t/h and for uncorrected Gaussian white noise 
we have ^ — > 6^/^^. By balancing the noise with the left-hand side, we find a; = 1; thus 
the driving field is marginal (?/ = 0), and the demagnetizing factor is relevant with scaling 
dimension w = 1. The demagnetizing factor k sets the characteristic maximum of the 
avalanche size and duration. The two-point probability function in the case c = can be 



found exactly 



25l | and the functional form of the resulting shape is, 



^V.^oAt\T)) = ^J^, ^, (10) 



and hence the scaling function is 



K, 0) = y . (11) 

Intuitively, the form of the scaling function leads to a flattening of the average shape at long 
times because k acts to cap the velocities. As it can be seen in Fig. [S3| this behavior is 
verified by our experiments. 

The effect of the driving field c, a marginal field at mean-field, can be calculated exactly 



at = using a result from Ref. 26|]: for an absorbing boundary condition at ^ = 0, 
the two-points probability function for the Brownian motion in a logarithmic potential is 
G'c,fc=o(V",t;e,0) = j^(4t)-(i+'^)e2^re-^'/(^*)(l + ©(e^)), where z/ = (1 - c)/2. Using this 
expression, the scaling form is reduced in magnitude but remains parabolic, with V(A, 0, C) = 
(1/8)(1 + C)A(1-A). 

In conclusion, we have used multivariable universal scaling functions to study average 
avalanche temporal shapes in Barkhausen noise. By analyzing thin films, we have bypassed 
the eddy current effects that have long frustrated complete comparison between theory and 
experiment. We unify the two rival mean field theories, calculate both the mean-field scaling 
function for the temporal average shape and the effects of the demagnetizing factor and the 
field rate. By utilizing optimal Wiener filtering techniques we allow an unbiased capturing 
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of the shapes, allowing us to report almost symmetric shapes, yielding excellent agreement 
with our theoretical predictions. 



Methods 

Sample preparation and experimental measurements 

A 1 yum thick ferromagnetic film with nominal composition NigiFcig (permalloy) is de- 
posited by magnetron sputtering from a commercial target, on a glass substrate covered 
by a 2 nm thick Ta buffer layer. The deposition is performed with the substrate moving 
at constant speed through the plasma in order to improve the film uniformity, with the 
following parameters: base vacuum of 10^'' Torr, deposition pressure of 5.2 mTorr (99.99 % 
pure Ar at 20 SCCM constant flow), and 65 W RF power. The deposition rate of 0.28 nm/s 
is calibrated with x-ray diffraction, which also confirms the polycrystalline character of the 
film. Quasi-static magnetization curves obtained with a vibrating sample magnetometer in- 
dicate isotropic in-plane magnetic properties with an out-of-plane anisotropy contribution, a 
behavior related to the stress stored in the film, and to the columnar microstructure 
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29|. 



Barkhausen noise time series are obtained using the inductive technique in an open mag- 
netic circuit. The sample has dimensions 10 mm x 4 mm. Sample and pickup coils are 
inserted in a long solenoid with compensation for the effects of the border. The sample is 
driven by a triangular driving magnetic field with amplitude high enough to saturate the film 
(± 15 kA/m). The driving field frequency is varied in the range 0.05 Hz-0.4 Hz. Barkhausen 
noise is detected by a 400 turn sensing coil (3.5 mm long, 4.5 mm wide, 1.25 MHz resonance 
frequency) wound around the central part of sample. A second pickup coil, with the same 
cross section and number of turns, is adapted in order to compensate the signal induced by 
the magnetizing field. The Barkhausen signal is then amplified and filtered with a 100 kHz 
low-pass filter, and finally acquired data at 0.2 and 4 MSamples/s. The data at different 
rates are statistically very similar after Wiener filtering (spurious peaks are present in the 
high frequency regime of the spectrum, independently of the rate), with very small differ- 
ences at short durations and small voltages. Wiener filtering is practically more efficient at 
large rates, for uncorrelated noise, but the noise amplitude remains larger than at low rates. 
Thus, the optimal use of this filter depends on the observables studied: For the distributions 
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in Fig. [S2l in order to gain more accuracy at short durations and small avalanches sizes by 
reducing the noise level, the low rate is used; For the average temporal shapes, in Fig. [SH 
being more important to capture variations of the signal at small timescales by a more accu- 
rate Wiener filter, the large rate is used. The time series is acquired just around the central 
part of the hysteresis loop near the coercive field, where the domain wall motion is the main 
magnetization mechanism 



23 



30| • By usin g a thin film, we have removed the confounding 



effects of eddy currents, whose time scale 



12| for the present sample is estimated to be 



Tp ~ 0.04 /is, much smaller than the avalanche durations studied. 



Data analysis 

To address the low signal from the thin film, we analyze the data by using Wiener 
deconvolution {l^. The form of the output signal is assumed to be of the form Vout{t) = {h-k 
(y + n)){t) where h(t) denotes the impulse response function, V(t) the original microscopic 
signal, and n{t) the background noise. Given an estimated impulse response Fourier series 
h{f), an estimated deconvolved noise power spectrum and a theoretically expected 

frequency spectrum for the deconvolved signal the filtered data V{t) is the inverse 

Fourier transform of 



V{f) = 



VoutU) 




kf) 


n{f)? + 





(12) 



Here, the impulse function is estimated by measuring the impulse response of the individual 
instruments composing the apparatus. We note that if the signal has a high resolution, a 
detailed knowledge of the impulse response function is needed to remove the effect of the 
instrument. As a signal function, we use a power fit ~ (see also Supplementary 

Information). To estimate the noise contribution |n(/)p, we measure the power spectrum 
of the instrumental noise, recorded without the material, maximizing over several runs. The 
Wiener deconvolution, shown in Fig. [ST] and in the inset of Fig [S2b . smoothens the signal, 
and removes spurious high frequency oscillations, due to the amplifier and filters used in the 
experiments. Importantly, this procedure also allows us to avoid the use of thresholds for 
defining the temporal extent of the avalanches, improving our estimates for both the scaling 
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exponents and the average avalanche shapes. 
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Figure Legends 




0.02525 0.0253 0.02535 

t(s) 

FIG. SI: Extracting average shapes from noisy data by Wiener deconvolution. Time- 
series data (jagged line) is traditionally separated into avalanches using a threshold Vth set above 
the instrumental noise (dotted blue line) - here breaking one avalanche into a few pieces. We 
instead do an optimal Wiener deconvolution (smoothed red curve, see text), allowing the use of a 
zero threshold (solid black line) which avoids distortions of the average shape and also gives more 
decades of size and duration scaling. Averaging over all avalanches with this duration gives us 
{V{t,T)) (dashed green curve). 
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lo' Duration T {\xs) lo' 2e-8 Voltage V (V) i-2e-7 

FIG. S2: Multivariable scaling; experiment vs. mean-field theory. In each main graph, the symbols represent data 
taken at various sweep frequencies Q for the external field. The lines represent mean-field predictions extracted from ABBM 
simulations (cf. Eq. [Gj with parameters c and k shown in the insets of (a) and (b), taking similar number of samples with 
the experiment; the sweep rate c oc f2 as predicted, and the weak dependence of the demagnetizing factor k can plausibly be 
attributed to changing experimental conditions. The distribution of avalanche sizes and durations shown in (a) and (b) depend 
on three variables (demagnetizing factor, sweep frequency, and S or T); no analytical form is available, but the agreement with 
the mean-field simulations is excellent. Avalanches of longer duration merge when ramped at a finite rate, leading to unusual 
rate-dependent critical exponents in mean-field theory; the smaller avalanches have P ~ S~'^ with r = 3/2 — c/2, while the 
briefer ones have P ~ with a = 2 — c. (Vertical scales are shifted for clarity; vertical units are arbitrary.) The limiting 

values of c = and c = 1 (corresponding to Q = 0.3) are shown. The mean avalanche size versus duration is shown in (c); for 
avalanches below the demagetization cutoff (S) ~ T^, while above the demagnetization cutoff (S) ~ T is observed, as expected. 
Power-law behavior of avalanche sizes versus durations has been observed before in experiments on systems with short-range 
interactions, but not before on materials believed to belong to the long-range or mean-field universality class (presumably 
because of eddy currents) . The observed systematic vertical shift with rate of the experimental curves is due to the presence of 
instrumental noise, not included in the simulations. The inset shows the power spectrum for f2 = 0.05 Hz, which has a behavior 
consistent with as predicted by mean-field theory. Finally, in (d) the the voltage probability distribution is consistent 

with the mean-field prediction P{V\k,c) oc V"^~^ exp (— fcy), another three-variable scaling form. The observed deviations at 
small voltages are due to the presence of instrumental noise, being approximately gaussian with tr « l8nV. 



Supplementary Information 



Beyond power laws: Universality in the average avalanche shape 

Stefanos Papanikolaou, Felipe Bohn, Rubem Luis Sommer, Gianfranco Durin, 
Stefano Zapperi, and James P. Sethna 



1 





- 38.5ns 

- 54.5ns 

- 77ns 

- 109ns 

- 154ns 
218.5ns 




(a) 



(b) 




0.2 0.4 t/T 0.6 Q.i 

(c) 




0.4 t/T 0.6 

(d) 



FIG. S3: Experiment and theory: Average shapes and scahng. (a) The temporal average 
avalanche shape for different avalanche durations T, rescaled to unit height and duration. The 
218.5 us curve, for example, represents the average over 658 avalanches, with durations between 
181 /US and 256 fis, binned in 50 equal time intervals, each of size 4.37 fis . Notice the symmetric, 
parabolic shape for short durations, and the flattening for longer durations. The small residual 
asymmetry remaining is plausibly attributable to amplifier-induced distortions that were not possi- 
ble to correct for. The impulse response functions of the individual instruments (amplifier, low-pass 
filter and sensing coils) were measured subsequent to taking the data, and used to deconvolve the 
resulting signal (cf. Supplementary Information), (b) Mean-field calculation for the average shape 
including a demagnetizing factor k as in Eq. [5j The shape is an inverted parabola, 4t/T(l — t/T) 
for short durations and fiattens at long durations. Here, k is 0.05 /us~^, slightly different from the 
value used in Fig. IS2[ This difference is expected since our analytical form is only available for 
c = 0; simulations at c, k given in the insets of Fig. [S2k and b respectively, yields good, albeit noisy, 
agreement with experiment. (c,d) Deviations from the inverted parabola. Note the quantitative 
agreement between experiment and theory. 

In this supplementary information, we 



(A) Explain in detail how we implemented the optimal Wiener filtering to extract the mag- 
netic avalanche sizes, durations, and temporal shapes from a noisy Barkhausen signal. 
We show that our methods not only allow for extraction of avalanche shapes, but provide 
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more faithful extraction of avalanche sizes, durations, and velocities from the signal, by 
comparing filtered and unfiltered analyses both of experiments and of simulations with 
added noise. 

(B) Further characterize our films, and discuss issues involved in the 3D-2D crossover, 
expected in the study of thinner films. 

(C) Explain the theoretical mapping between the two discussed mean-field theories and 
describe the approximation made in the previous average shape calculation (which 
gave a sine scaling function rather than an inverted parabola) . In order to demonstrate 
the technique, we show in detail how to calculate the average shape ior k — c — 0. 

A. Optimal Wiener filtering for extracting Barkhausen avalanches 

In this publication we focus on thin magnetic films to avoid the complexities of eddy cur- 
rents found in ribbons. Because of the reduced thickness, the induced signal is significantly 
smaller, substantially decreasing the signal-to-noise ratio. The larger number of turns of 
pick-up coil compensates for the reduced voltage, but also introduces spurious background 
effects, mainly in the high frequency regime of the spectrum. Such effects clearly demand a 
more sophisticated extraction of the signal from the noise. We found that the most effective 
results are obtained by applying an optimal Wiener filter, using independent measurements 
of the spectral properties of the background noise. 

le-08 



le-10 
le-12 
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FIG. S4: The average power spectrum of the signal. The power spectrum at high frequencies, 
averaged over 200 runs, is dominated by background noise effects which alter the scaling behavior 
at intermediate sizes and durations. 
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(a) (b) 

FIG. S5: The power spectrum of the instrumental noise after one sweep and the im- 
pulse response. The black points in (a) are the power spectrum of one Barkhausen run, out of 
200 total used in this paper, deconvolved with the function shown in (b), based on the measured 
impulse response function for the instrument. The red points in (a) denote the background decon- 
volved noise spectrum, which (after deconvolution) resembles white noise at low frequencies. The 
instrumental response in (b) is composed of a low-pass filter response, which is a pure exponential 
decay and a sensing coil response, which amounts to an oscillatory exponential decay. 
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FIG. S6: The effect of instrumental noise on the average shapes. The effect of the back- 
ground noise oscillations on the form of the average shapes is easily observed in the individual 
avalanches (Fig. 1 in the main text), and dramatic in the average avalanche shapes for fixed (non- 
binned) durations (above). After filtering, the oscillations are reduced significantly. 

Optimal Wiener filtering is performed in Fourier space, using the available information 
about the spectral properties of the signal and the background noise. The signal is pre- 
dominant in the low-frequency portion of the spectrum where, in general, it shows a power 
law behavior at intermediate frequencies, going to a constant at low frequencies. Since we 
are interested in filtering the high-frequency regime of the spectrum, we are very careful 
about estimating the expected signal S{u) at the intermediate frequencies, where scaling 
takes place. For this purpose, we use a power law fitting function S{ijj) ~ for our 
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FIG. S7: The power spectrum and Wiener deconvolution. (a) The black squares denote 
the power spectrum of 200 Barkhausen runs for the rate 0.05Hz, deconvolved with the impulse 
response function of Fig. ISSf b). The red curve is a power law fit to the scaling regime of the 
averaged signal, used in the Wiener deconvolution method. The purple circles denote the average 
of the background noise spectrum, deconvolved in the same way as the signal, over 200 independent 
measurements. Finally, the blue curve denotes the estimated noise spectrum, used in the Wiener 
deconvolution method, which is composed of a white-noise component located at the average 
amplitude of the background spectrum in the low-frequency regime, and the spurious isolated 
peaks at high frequencies. Since the positions of the spurious peaks move from run to run, we 
include them in the noise-estimate whenever the background-noise spectrum of a single run has 
larger than 3(T„, where (T„ is the variance of the noise. The form of the signal is plagued with 
high frequency features that come from artificial background noise elements, (b) After the Wiener 
deconvolution, the resulting filtered signal's spectrum averaged over all runs (green diamonds) 
reduces significantly the distorted features, recovering a large part of the scaling behavior, which 
is consistent with the mean-field prediction (^(a;) ~ w^^). 



expected signal (cf. Fig. [SW a)). In this way, the very low frequency regime remains prac- 
tically unfiltered (since the spectrum intensity of the signal is much larger than that of 
the noise (cf. Fig. ISST a). ISW a)). The background noise, after deconvolving the instrumen- 
tal response (cf. Fig. ISSt^ a). [SS^ b)). is white at low frequencies (cf. Fig. ISSt^ a) and ISW a)). 
The high-frequency portion of the signal and noise are affected by nonlinear amplifier arti- 
facts (cf. Fig. ISST b)). which produce a series of resonance peaks (cf. Fig. [SH ISST a)). Their 
effects on the time signal are spurious oscillations which persist also in the average avalanche 
shape, as seen in Fig. [S6l We use this information to set our estimate for the noise spectrum 
as a sum of two components: i) a white-noise component, being the mean of the Gaussian 
noise (~ lO^^^V^Hz^^), shown in Fig. ISW a). ii) Spurious peaks in the spectrum, with am- 
plitudes larger than three times the white noise's variance, that are observed to contaminate 
the signal at high-frequencies. Since these noise peaks are quite sharp and slightly shift from 
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run to run, we use the maximum of these 3a events over our 200 independent magnetic field 
sweeps as our estimated noise spectrum (blue curve in Fig. [STFa)). 

With this information on the instrumental noise and on the signal, we apply Eq. 12 of 
the main text to the FFT of the time signal (cf. ISTlfa)). and perform the inverse FFT to get 
the filtered time signal (cf. [HW b)). The average avalanche shapes now appear smooth (see 
Fig. 2 in the main text). 

We test our optimal Wiener filtering by studying its effects on avalanche sizes, durations, 
and voltages in Fig. [Si We find that the method appears to have significant advantages over 
both analyzing the unfiltered data and the alternative method of thresholding for extracting 
smaller avalanches from data (Fig. [581(a) and (c))). We also compare to simulated data 
(cf. Eq. 6 of the main text) - we add Gaussian noise to match that of the experiment after 
deconvolution, and then we convolve with the experimentally measured impulse function, 
to test for potential systematic errors due to the Wiener filtering. We find that the smaller 
and shorter avalanches and the low voltages (below the noise) are indeed distorted by the 
process, but that the recovery of large avalanches is significantly enhanced by the Wiener 
filtering process. 

We also use simulated data to explore possible systematic effects on the average avalanche 
shapes. After adding the noise and convolving with the impulse response function of 
Fig. ISST b). the average shape develops a significant asymmetry that resembles the one seen 
in the experiments, shown in Fig. [S9l After filtering, this asymmetry disappears completely, 
since in this case we know the form of the impulse response exactly. In the case of the actual 
experimental shapes (cf. Fig. 2 of the main text), even though filtering reduces the present 
small asymmetry, it does not disappear completely, presumably because the true impulse 
function involves non-linear effects that are hard to identify and correct for. Also, note that 
the noisy signal includes a systematic correction near the edges of the shape (t — > and 
t — )■ T), due to the presence of background noise. 

B. Films and dimensionality 

Here we characterize our experimental film in more detail, and explain why it is not 
surprising that (as clearly shown in Figs 2 and 3 in the main text) they exhibit critical 
behavior consistent with the mean-field theory describing three-dimensional magnets. 
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FIG. S8l Optimal Wiener filtering effects, tested on experimental data and simulations with added noise. For 
each pair, the left panel is experiment and the right panel is simulated ABBM mean-field data with k = 0.08, c = 0.15, and 
added noise roughly matching that of the experiment (at rate 0.05 Hz), after deconvolution. (a) and (b) show the avalanche 
size distributions; (c) and (d) show avalanche duration distributions, and (e) and (f) show voltage distributions. For the 
experimental size and duration distributions (a) and (c), the optimal Wiener filtering (green) appears to allow the extraction of 
almost one decade more than either the raw data (red) or data analyzed with the noise threshold depicted in Fig. 1 of the main 
text. The theoretical size and duration distributions (b) and (d) show the true ABBM distributions (black), the distributions 
analyzed with added noise (red), and the optimally filtered simulations (green). Here we see that the filtering allows us to 
recover many large avalanches that split apart under the noise, but that the distribution of small and short avalanches is 
notably distorted. The voltage distributions (e) and (f) are distorted by the noise and filtering only at low voltages. 

Avalanches whose spatial extent is large compared to the thickness of the film should 
eventually be described by a two-dimensional theory. Indeed, for very thin films (below 



200 nm), new universality classes do appear [l-|3|. For relatively thick polycrystalline rib- 
bons (down to 45/im), two of us have previously shown that Barkhausen noise shows 
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FIG. S9: Simulated avalanche shapes after Wiener deconvolution for the ABBM model 
before and after filtering. Two durations of 109 fis and 218 //s shown (where each simulation 
time-step is taken to be 0.25 /is as in the experiments), before and after filtering the ABBM time 
series with added Gaussian noise and convolved with the experimentally observed filters, (a) The 
differences between the ABBM and noisy shapes is significant, but after filtering the behavior is 
much improved, (b) After adding the noise and convolving with the impulse response function of 
Fig. ISSf b) , the average shape develops a significant asymmetry that resembles the one seen in the 
experiments. After filtering, this asymmetry disappears completely, since in this case we know the 
form of the impulse response exactly. In the case of the actual experimental shapes (cf. Fig. 2 of 
the main text), even though filtering reduces the present small asymmetry, it does not disappear 
completely, presumably because the true impulse function involves non-linear effects that are hard 
to identify and correct for. Also, note that the noisy signal includes a systematic correction near 
the edges of the shape (t — )• and t — )• T), due to the presence of background noise. 



the mean-fiejd behavior expected for 3D interface depinning problems with strong dipolar 
interactions [5] - the observed avalanches are small enough compared to the thickness, thus 
they should not be affected by the asymptotic crossover to 2D behavior. Our sample has 
an intermediate thickness of l/im, and its polycrystalline state is confirmed by the x-ray 
structure shown in Fig. ISlOi The magnetization curves of Fig. ISllI show the typical shape 
exhibited by thick films with a weak out-of-plane anisotropy, due to the stress stored in the 
film during the deposition and to the columnar microstructure Q]. The in-plane magnetic 
properties are essentially isotropic, with domains of dense stripes having widths of the order 
of the thickness, as confirmed by performing MFM measurements. Such behavior is known 
to be absent in the experiments showing 2D behavior iMSj. 
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FIG. SIO: High angle x-ray diffraction pattern confirming the polycrystalline state of 
the film. The (111) and (200) permalloy peaks are identified at 26 ~ 44.2°, and 29 ~ 51.5°, 
respectively. 
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FIG. Sll: In-plane normalized magnetization curves measured by VSM, obtained along 
two different perpendicular in-plane directions = 0° and (j) = 90°. The angle (j) is 
between the applied magnetic field and direction defined by the motion of the substrate during 
the deposition. The magnetization curves indicate isotropic in-plane magnetic properties with an 
out-of-plane anisotropy contribution, a behavior related to the stress stored in the film during the 
deposition, and to the columnar microstructure 



Connections betw^een the ABBM & Shell models 



Here we derive in more detail the connection between the two mean-field theories (the 
ABBM and shell models), and then implement in detail the calculation of the paraboloic 
average avalanche shape for the case A; = c = 0. 



The equation for the ABBM model (Eq. 2 o 
simpler equation after taking a time derivative j?!: 



the main text) can be transformed to a 



r§^^...(x). 



(13) 
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The noise term has delta-correlations in magnetization 



{w{X)w{X'))^ = 2D5{X - X'). (14) 

Also, dX/dt = V and = |. 

Eq. [13] belongs to a class of equations of the form 

^ = f{V)+aX), (15) 

with 

mnx')) = cTis{x-x'). (16) 

The units of [^] = 1/t, so the units of [ax] = y/X/t. We interpret Eq. [T5] as a rule for 
incrementing V and t in discrete steps of size AX = e: 



Vx+e = Vx + ef{V) + E{X), (17) 







where now H(X) is an uncorrelated Gaussian variable of standard deviation as = \feax 
This can be seen explicitly; cr= represents the RMS width of a random walk with noise S of 
duration e: 

+ e) - x(X))2) = I ' ci/i'(e(/x)e(/x')) = ^4, (18) 

using Eq. [T6l 

Now let us consider finding V^(t + At), where At is large compared to tjV (so many 
m-steps are averaged over) but small enough that V doesn't change significantly during the 
interval [8]. We do this by repeating Eq. [TTl N = V ^t/e times. Then 

X{t)+V At 

V{t + M)^V{t) + f{y)V At+ S(X). (19) 

X=X{t) 

The sum has X terms, so is a Gaussian uncorrelated variable with standard deviation given 
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by 



y/Nas = \l axV^ = VVAtax- (20) 



If we define a new discrete Gaussian variable Z of standard deviation az = VAtax, then 
multiplying Z by y/V mimics the sum, so we find 

V{t + At) = V{t) + Vf{V)At + Z{t)Vv. (21) 



If we now change to continuous time, then (analogous to Eq. [T8l) az = a^v At, so 

dV 



VfiV) + Vvat), (22) 



yielding Bertotti's result quoted in the main text. Here, ((t) satisfies 



{amt')) = ^^S{t-t') = aj,S{t-t'). (23) 



From Eq. EH the units of [C] = y/X/t^; from Eq. [23l the units of [(CC)] = X/t^, matching 
the units of [(Tx?[S{t)] = {X/e)[l/t]. 

A nice physical idea exists behind Eq. [221 Namely, the motion dX/ dt is always forward, so 
it always uncovers "new" random forces as the front progresses. Thus the noise as a function 
of time must be delta-function correlated (the new sites being passed are uncorrelated with 
sites passed at earlier times). The key question then becomes what the strength of the delta 
function is: how does the noise felt by a particle moving fast over a random environment 
differ from that felt by a particle moving slowly? Basically, one is summing over a number 
of random variables proportional to V , so the noise is increased by a factor of square root 
of this number, \/V . 

All the above manipulations have been performed using the convenient Stratonovich 
interpretation. In those cases where we have calculated average shapes using the more 
physical Ito calculus, the results have remained unchanged, up to dimensionless constant 
factors. We conjecture that the two methods, although different in detail, do not differ in 
their predictions for universal quantities. 
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The temporal average avalanche shape as a first return of a random walk to the origin 



Here, we show explicitly, in the simplest possible example, how one can use the concept 
of the first return of a random walk to the origin in order to calculate the form of the average 
shape's scaling function. Let's consider the (no-bias, constant-field) hmit where c — > and 
A; ^ 0, where it becomes equivalent to the so-called shell model: 

^ = ^«*)- P*) 

as described in the main text. Changing variables to x — VV, this equation becomes an 
ordinary random walk. We solve for an initial condition at e with absorbing boundary 
conditions at the origin, which yields for e — > 0, 

Grw{x, t; e, 0) = ^^'/'e xe-^'/^''\ (25) 
after using the method of images. In the original V coordinates it becomes: 



G.(y,t;e,0) = ^«"'')y-!^^-°) . (26) 

The average avalanche shape is defined as the average of the signal V{t,T) at time t 
for an avalanche of duration T. Given the probability distribution Gs, the probability 
of a signal V{t,T) at time t is Gs{V,t;e,0)Gs{e,T;V,t) = GsiV,t;e,0)GsiV,t;e,T) = 
Gs{V,t; e,0)Gs{V,T — t; e,0), using time-reversal and time-translation invariance of our 
random walk, where Gs has an absorbing boundary at = e. In other words, the average 
of V{t, T) naturally is 

{V{t,T)) = dVVGs{VAe,Q)Gs{V,T -t-e^Q) 



e^o f^°^ dVGs{V,t;e,0)GsiV,T -t;e,0) 
f 



'^T{Ul'^]]+0{e), (27) 
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and therefore, the shape is an inverted parabola. 
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